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Abstract. In the framework of percolation of strings, the transverse momentum distributions in AA and 
hh collisions at all centralities and energies follow a universal behavior. The width of these distributions 
is related to the width of the distribution of the size of the clusters formed from the overlapping of 
the produced strings. The difference between the distributions for baryons and mesons originates in the 
fragmentation of clusters of several strings which enhance the particles with higher number of constituents. 
The results agree with SPS and RHIC data. The predictions for LHC show differences for baryons compared 
with RHIC. At LHC energies we obtain also a high pt suppression for pp high multiplicity events compared 
with pp minimum bias. 

PACS. 25.75.Nq - 12.38. Mh - 24.85+p 



Multiparticle production can be described in terms of 
color strings stretched between the partons of the projec- 
tile and target. These strings decay into new ones by sea 
q — q pair production and subsequently hadronize to pro- 
duce the observed hadrons. The color in these strings is 
confined to a small area in transverse space: Si = nr^ with 
To ~ 0.2 — 0.3 fm. With increasing energy and/or atomic 
number of the colliding particles, the number of exchanged 
strings grows and they start to overlap , forming clusters, 
very much like disks in two-dimensional percolation the- 
ory. At a certain critical density, a macroscopical cluster 
appears, which marks the percolation phase transition [TJ 
2 . For nuclear collisions, this density corresponds to the 
value of r\ = Ns-§j of r\c — 1-18 — 1.5 (depending of the 
type of the employed profile functions) where Ns is the 
number of strings and Sa corresponds to the overlapping 
area of the nuclei. A cluster of n strings behaves as a sin- 
gle string with and energy-momentum that corresponds 
to the sum of the energy-momentum of the overlapping 
strings and with a higher color field, corresponding to the 
vectorial sum of the color fields of the individual strings. In 
this way, the multiplicity < fi n > and the mean transverse 
momentum squared < p\ >„ of the particles produced by 
a cluster are given by 



EqsQ] transform into analytical ones [5] in the limit of 
random distributions of strings 



VnS n 2 nS\ 2 

< fi >i <p T >„= W — < p T >i 

. w 

where < /x >i and < p\ >i stand for the mean multiplic- 
ity and mean p\ of particles produced in a single string. 



< n >= NgF(r)) < m >i 



<p z T >-- 



< Pt >i 
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where F(n) 



If we are interested in a deter- 



mined kind of particle i, we will use < /i >u, < p\ >u 
< [i > n i and < p\ > n i for the corresponding quantities. 
The transverse momentum distributions can be written as 
a superposition of the transverse momentum distributions 
of each cluster, g{x,pr), weighted with the distribution 
of the different tension of the clusters, i.e the distribu- 
tion of the size of the clusters, W(x) [3 : ,.,4j . For g(x,px) 
we assume the Schwinger formula g{x,pr) = exp{—p\x) 
and for the weight function W(x), the gamma distribution 
W{x) — Yjj?j(~fx) k ~ 1 exp(—jx) where 



< x > 



< or > - < x 
< x > 2 



x is proportional to the inverse of the tension of each clus- 



k is propor- 



ter, precisely x = 1/ < p\ > T 



nSi <?4> ' 

tional to the inverse of the width of the distribution on x 
and depends on 77, the density of strings. Therefore, the 
transverse momentum distribution f(pT,y) is 



dN 
dp%,dy 



f(pr,y) 



dxW(x)g(pT, x) 
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dN k - 1 



dy k < p r >i 



■F(v)- 



(1 



k<p'i,>u 



(3) 



Eq.Q is valid for all densities and type of collisions. It 
only depends on the parameters < pt >h and k. At low 
density, there is no overlapping between strings and there 
are no fluctuations on the string tension , all the clus- 
ters have one string. Therefore, k goes to infinity, and 

2 

f(pr,y) — exp(- <p p / >i ). At very high density rj, there 

is only one cluster formed by all the produced strings. 
Again, there are no fluctuations, k tends to infinity and 



f(j>T,y) - exp(- 



<K>i 



) . In between these two limits, 



k has a minimum for intermediate densities correspond- 
ing to the maximum of the cluster size fluctuations. This 
behavior of k with rj [311] is related to the behavior with 
centrality of the transverse momentum [5] and multiplicity 
fluctuations [6]. We observe from eq.© that 



dlnf{p^,y) 
dlnpT 



2F( V ) 



Pt 



(1 



) < p T >U 
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As p r — > 0, eq.(H]) reduces to 2F(i))p T w jjji e f or l ar g er 

Pt it becomes —2k for all particle species. As < p T >i p >< 
Pt >ik>< Pt >i7r, the absolute value is larger for pions 
than for kaons and for protons, in agreement with exper- 
imental data [3j. 

The nuclear modification factor , defined as Raa(pt) = 
"' A /N co u ^ P ^ y , reduces to the following expression at 



dp^dy 

Pt 



(we use eq.2 for ^): 
Raa(0) ~ ( 



(5) 



where rj' and rj are the corresponding densities for nucleus- 
nucleus collisions and pp collisions respectively, rj' > r/, 
thus F(r]') < F(r]). As pt increases we have 



Raa{Pi 



Pt 



<p T >n 



1 + F(rf) 



FW )2 
F( V ) ) 



(0) 



<p t >h 

and Raa increases with pt up to a maximum value. At 



Raa 



2(k(r,) 

Pt 



-fc'(V)) 



(7) 



At high density k'{rj') > k(rf) and suppression of pt oc- 
curs. 

The universal formula © must be regarded as an an- 
alytical approximation to a process which consists on the 
formation of clusters of strings and their eventual decay 
via the Schwinger mechanism. We do not claim to have 
an alternative description valid at all pt- It is well known 
that jet quenching is the working mechanism responsible 
of the high pt suppression. This phenomena is not in- 
cluded in our formula which was obtained assuming a sin- 
gle exponential for the decay of a cluster without a power 



like tail. Our work must be considered as an interpolating 
way of joining smoothly the low and intermediate pt re- 
gion with the high pt region. The suppression of high pt 
has correspondence with a modification in the behavior 
at intermediate pt and this is what we study. Although it 
describes many of the observed features of experimental 
SPS and RHIC data, it is not able to explain the differ- 
ences between antibaryons(baryons) and mesons. In fact, 
the only different parameter between them in formula© 
is the mean transverse momentum of pions and protons 
produced by a single string < p T >iir and < p\ >\ v re- 
spectively. This only causes a shift on the maximum of 
Raa but keeps the same height at the maximum contrary 
to the observed. However, in the fragmentation of a cluster 
formed by several strings, the enhancement of the 

production of antibaryons(baryons) over mesons is not 
only due to a mass effect corresponding to a higher tension 
due to higher density of the cluster (factor F(t\)\vl front of 
p\ in formula © ). In fact, the color and flavor proper- 
ties of a cluster follow from the corresponding properties 
of their individual strings. A cluster composed of several 
quark-antiquarks (q — q) strings behaves like a (Q — Q) 
string, with a color Q and flavor composed of the flavor 
of the individual strings. As a result, we obtain clusters 
with higher color and differently flavored ends. For the 
fragmentation of a cluster we consider the creation of a 
pair of parton complexes QQ [7] . After the decay, the two 
new QQ strings are treated in the same manner and there- 
fore decay into more QQ strings until they come to objects 
with masses comparable to hadron masses which are iden- 
tified with observable hadrons by combining into them the 
produced flavor with statistical weights. In this way, the 
production of antibaryons(baryons) is enhanced with the 
number of strings of the cluster. As an example, in fig.l we 
show the results for the decay of a color octet cluster and 
color sextet formed by two 3 — 3 strings [7] , were there is 
a large enhancement of antibaryons(baryons). Notice that 
the enhancement of strangeness with zero baryon number 
is smaller [ll[8j[9|. We observe that the additional anti- 
quarks (quarks) required to form an (anti)baryon are pro- 
vided by the antiquarks (quarks) of the overlapping strings 
which form the cluster. In this way, the recombination and 
coalescence ideas |10j are naturally incorporated in our ap- 
proach. In order to take this into account in our formulas, 
we must modify the eqs.fl]),© an d ©• For (anti)baryons 
we will consider the multiplicity per unit of rapidity to be 
instead of the first equation of [2j 



(8) 



fitting the parameter a to reproduce the experimental de- 
pendence of the pt integrated p spectra with centrality 
[IT] . The result is shown in fig. 2 and the obtained values 
are a = 0.09 and — ^ = /j>ip and fiiw are the mean 
multiplicity of a single string for anitprotons and pions 
respectively. It is observed that when an antibaryon is 
triggered, the effective number of strings is N$ +a instead 
of Ns- This means that the density rj must be replaced 
by 1Tb ~ -^Ys 7 ?- The (anti)baryons probe a higher density 
than mesons for the same energy and type of collision. On 
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the other hand, from constituent counting rules [12j it is 
expected that the power-like px behavior for baryons is 
suppressed in one half more than mesons. Therefore, in 
((3J we must use for (anti)baryons tjb and the correspond- 
ing functions F{j]b) and ks = k{y\B) + \- Since tjb > V 
we have F(t]b) < F(rj) and k(r)B) > Mjj). For periph- 
eral collisions, Ng is smaller than for central collisions 
and r\B is more similar to r\. Therefore, the differences be- 
tween the transverse momentum distributions are smaller 
as it is shown by the experimental data. In order to com- 
pute the different 77s for different centralities we use the 
Monte-Carlo code of references 7 , 8 . Their values at RHIC 
and LHC and their corresponding k's are tabulated below. 
Note that the values of k come from the universal function 
which gives the shape of the dependence of k on 77 J3j[4] . 

In our approach, we should use a different a for baryons 
and antibaryons as far as the increase with centrality is 
slightly different for both, depending also on the specific 
kind of baryon(antibaryon). In order to obtain a single 
formula we considered strings of the same type. However, 
usually two types of strings are considered. Strings qq — q 
or q — qq which stretch a diquark of the projectile (target) 
with a quark of the target(projectile) and strings q — q or 
q — q linking quarks and antiquarks. The fragmentation 
of the strings of the type qq — q or q — qq favors the pro- 
duction of baryons over antibaryons in the fragmentation 
regions of the projectile and the target. Therefore our re- 
sults should be limited to the central rapidity region. On 
the other hand, we have determined a from the depen- 
dence of ^ with centrality. The large error data translate 
into uncertanties in a of the order of 20% which is of 
the same order of the differences in the parameter a for 
baryons and antibaryons. The uncertanties in a induce 
uncertanties in the determination of tjb and hence in ks] 
these uncertanties are however negligible (less than 5% 
even at the highest centrality). We conclude that our sin- 
gle formula using the same a for baryons and antibaryons 
is a good approximation. 

On the other hand, we have assumed an additional dif- 
ference between baryons and mesons, ks = k(rjB) + 1/2, 
from the high px behavior of eq.3. This could seem in- 
consistent with our approach that is limited to low and 
intermediate px- As we have said above, our goal is not 
to give a full description of data, including high pt but 
an alternative description of intermediate pt suppression. 
In this way, this factor can be seen as a boundary condi- 
tion imposed to our approach. Notice that this factor is 
independent of the centrality and therefore its influence 
on the ratios Rcp and Raa is not large. In fact, if we use 
for ks the values without the factor 1/2, Rcp is smaller 
at pt = 2 GeV/c and px = 10 GeV/c in a factor 1.13 and 
1.25 respectively at RHIC energies. This small reduction 
does not spoil the agreement with data although a better 
agreement at high px is obtained with the 1/2 factor. 

In fig 3., we show our results for the ratio Rcp in Au- 
Au collisions defined as usual, for (p +p)/2 (dashed line) 
and neutral pions (solid line), compared to the PHENIX 
[TT] experimental data. We find a good agreement. We 
also show the LHC prediction for pions and antiprotons. 



There is no change for pions but on the contrary, the differ- 
ence between antibaryons and pions is enhanced. In fig. 4 
we show the px dependence of the ratio for peripheral 
(dashed line) and central (solid line) Au-Au collisions to- 
gether with experimental data [TT]. LHC predictions are 
also shown. The difference only appears for central colli- 
sions. 

In fig. 5 we show our results for the modified nuclear 
factor Raa for pions (solid lines) and protons (dashed 
lines) for peripheral and central collisions together with 
the experimental data for pions(TT3"]). Again, the main 
difference arises for central collisions. In fig. 6, this differ- 
ence between pions and protons is compared at RHIC and 
LHC energies. The Cronin effect becomes larger at LHC 
for protons, contrary to some expect at ions 1 1 4j . 

The good agreement obtained with the experimental 
data can be understood as two combined effects:the larger 
string tension of the cluster and the formation of strong 
color fields [15] and the way of fragmentation of the clus- 
ters which enhances (anti)baryon over mesons similarly to 
recombination models [10!. Both effects are widely recog- 
nized as working physical mechanisms at high densities . 
Both effects are naturally incorporated in the percolation 
of strings approach. 

The shape of Raa and Rcp has nothing to do with the 
nucleon structure of the nucleus an it depends essentially 
on the string density. One can wonder whether in pp at 
LHC energies can be reached enough string densities to 
get a high px suppression. In fig. 7 we answer this ques- 
tion. It is plotted the ratio Rcp between the inclusive pp 
going to tt, k and p cross section for events with a mul- 
tiplicity twice higher than the mean multiplicity and the 
minimum bias cross section. It is observed a suppression 
for pt larger than 3 GeV/c. However, the reached string 
density in pp collisions will not be enough to suppress the 
back to back jet correlations as observed at RHIC energies 
for Au-Au central collisions. A straightforward evaluation 
following [3], gives a maximum of absorption of Apx = 2 
GeV/c for a number of strings of the order of 25-30 corre- 
sponding to events with a multiplicity twice the minimum 
bias multiplicity. 

We thank Ministerio de Education y Ciencia of Spain 
under project FPA2005-01963 and Conselleria de Edu- 
cation da Xunta de Galicia for financial support. We thank 
N.Armesto, C.Salgado and Y.Shabelski for discussions. 



References 

1. N. Armesto, M. A. Braun, E. G. Ferreiro and C. Pajares, 
Phys. Rev. Lett. 77 (1996) 3736; M. Nardi and H. Satz, 
Phys. Lett. B 442 (1998) 14. 

2. M. A. Braun and C. Pajares, Phys. Rev. Lett. 85 (2000) 
4864; Eur.Phys.J C85,349 (2000)'. 

3. J. Dias de Deus, E. G. Ferreiro, C. Pajares and R. Ugoc- 
cioni, Eur. Phys. J. C 40 (2005) 229. 

4. C. Pajares, Eur. Phys. J. C 43 (2005) 9; J. Dias de Deus 
and R. Ugoccioni, Eur. Phys. J. C 43 (2005) 249. 

5. K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 
757 (2005) 184; E. G. Ferreiro, F. del Moral and C. Pa- 



4 



L. Cunqueiro 1 et aL: Title Suppressed Due to Excessive Length 



Centrality | r] | k | 773 | fcs 
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3.97 


5.17 


4.82 


Au-Au 60 - 92% 


0.9 


3.58 


1.24 


4.16 
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0.6 


3.55 


0.75 


4.06 


PP 


0.4 
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0.48 


4.07 
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PP 


0.72 
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0.92 


3.90 



Table 1. Density of strings and the corresponding k values for 
mesons and barions at RHIC and LHC. 
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Relative to ;i~mean particle numbers 




Fig. 1. Hadronic content of the decay of octet and sextet 
strings relative to the number of n~ . 
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Fig. 4. Antiproton to neutral pion ratio as a function of pr Fig. 6. Nuclear modification factor for 0-10% central pions 
for for 0-10% (solid) and 60-92% (dashed) centrality bins com- (solid) and protons (dashed) at RHIC (black) and LHC (blue) 
pared to PHENIX data. LHC predictions in blue. 




qt 



Fig. 5. Nuclear modification factor for neutral pions (solid) 
and protons (dashed) for 0-10% central and 80-92% peripheral 
bins compared to PHENIX data. 



Fig. 7. Central to peripheral ratio for pp collisions at LHC. 
Black: ir , Red: kaons, Blue: p 



